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We study three-body nonleptonic decays B — > VVP by introducing two-meson distribution ampli- 
tudes for the vector-pseudoscalar pair, such that the analysis is simplified into the one for two-body 
decays. The twist-2 and twist-3 (j>K two-meson distribution amplitudes, associated with longitudi- 
nally and transversely polarized <j> mesons, are constrained by the experimental data of the r — > (f>Kv 
and B — > (pK'y branching ratios. We then predict the B — * (pK'y and B — + 4>4>K decay spectra in the 
<f>K invariant mass. Since the resonant contribution in the (j>K channel is negligible, the above decay 
spectra provide a clean test for the application of two-meson distribution amplitudes to three-body 
■ B meson decays. 
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Viewing the experimental progress on three-body nonleptonic B meson decays 0, 0, it is urgent to construct 
a corresponding framework. In j3( we have proposed a formalism based on the collinear factorization theorem in 
perturbative QCD (PQCD), in which new nonperturbative inputs, the two-meson distribution amplitudes, were 
introduced 3 . On one hand, a direct evaluation of hard kernels for three-body decays, which contain two virtual gluons 
at lowest order, is not practical due to the enormous number of diagrams. On the other hand, the region with the two 
| gluons being hard simultaneously is power-suppressed and not important. Therefore, the new nonperturbative inputs 
^\ . are necessary for catching dominant contributions in a simple manner. In our formalism the collinear factorization 
' formula for a B — > h\h-2h% decay amplitude is written, in general, as 

-xt- : 

O. M = $s ® # <8 <8 $fc s , (1) 

"xl" ■ 

where &B,h 3 are the B,h^ meson distribution amplitudes, <&hih 2 the \i\h2 two-meson distribution amplitude, and 
■ ® represents the convolution in longitudinal momentum fractions x. <&hih 2 and ^i la include not only the twist-2 
Oh! (leading- twist), but two-parton twist-3 (next-to-leading-twist) components. The computation of the hard kernel H, 
basically the same as in two-body B meson decays, is restricted to leading order in the coupling constant a s so far. 
I There are two types of factorization theorems: collinear factorization |a,SSS0 an d kx factorization [Tol [TT| . 
For a comparison of the two types of theorems, refer to Collinear factorization works, if it does not develop 

an end-point singularity from x — > 0. If it does, collinear factorization breaks down, and kx factorization is more 
. , appropriate. It has been known that collinear factorization of charmed and charmless two-body B meson decays 
suffers the end-point singularities fT3 . T his is the motivation to develop the PQCD formalism for two-body B meson 
decays based on kx factorization |!5lll6llT7| . This approach has been shown to be infrared-finite, gauge- invariant, and 
consistent with the factorization assumption in the heavy-quark limit |l8t Il9l l20|. For three-body B meson decays, 
the end-point singularities are smeared by the two-meson invariant mass [3j, and collinear factorization in Eq. 
holds. Moreover, it has been demonstrated that both nonresonant contributions and resonant contributions through 
two-body channels can be included by means of an appropriate parametrization of &h 1 h 2 13 ■ 

One of the challenges in the studies of three-body heavy meson decays is the evaluation of the matrix elements 
for heavy meson transition into two hadrons. There are already several theoretical approaches to this subject in the 
literature. The naive factorization [2l| for three-body B meson decays has been adopted in , in which the B meson 
transition into two hadrons was simply parameterized by a power-law behavior and then fit to experimental data. 
The matrix elements for the above transition were calculated using the pole model |23l |2J, |25j , in which intermediate- 
state decays into two hadrons were described by effective weak and strong Lagrangians. The naive factorization has 
been improved in a so-called QCD-factorization framework |26|. However, only the current-produced amplitudes, 
i.e., those which can be expressed as products of two form factors in the factorization limit, were studied. The 
challenging subject of the B meson transition into two hadrons was not addressed [2(j. Compared to the above 
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methods, our approach does not rely on the naive factorization, since the nonfactorizable contribution is taken into 
account through nonfactorizable hard kernels. It is complete in the sense that various topologies of amplitudes, such 
as the B meson transition into two hadrons and the current-induced one, are analyzed in the same framework. It is 
also more systematic, because sub-leading corrections can be evaluated order by order in a s and power by power in 
the ratios w/rriB and m^/m^, where w is the invariant mass of the two-meson system, and tub (w/^) the B (/13) 
meson mass. 

In [3| we have applied Eq. |QJ to the modes, in which both hi and hi are pseudoscalar mesons P. The modes 
with h\ being a vector meson V and hi a pseudoscalar meson P have been observed recently |27| . Hence, we shall 
extend our formalism to three-body decays involving the B — > VP transition, taking B — > <jxj>K as an example. 
We shall first define the 4>K two-meson distribution amplitudes, which are more complicated than the PP ones. A 
simple parametrization is then proposed, and constrained by the experimental data of the r — > (j)Kv and B — > <f>K^f 
branching ratios. Afterwards, we predict the decay spectra of the B — > (f>K"f and B — ► <jxj)K modes in the <pK invariant 
mass. The resonant contribution through the <j)K channel is expected to be negligible: the i^i(1650), ^(1770), and 
-ftT(1830) mesons decay into the 4>K pair with the branching ratios not yet available in 28]. Therefore, the above 
spectra provide a clean test for the application of two-meson distribution amplitudes to three-body B meson decays. 

Label the momenta of the <p and K mesons from the B meson transition as Pi and P2, respectively. The B meson 
momentum Pg and the total momentum of the <j)K pair, P = Pi + Pi, are chosen, in the light-cone coordinates, as 

P B = ^|(1,1,0 T ), P=^(l^, 0r ), (2) 

with the variable r\ = w 2 /m B . Define C = Pi I ' P + as the (f> meson momentum fraction and = m^/ms as the <fi 
meson- P meson mass ratio, in terms of which the other kinematic variables are expressed as 



P+ = (1 - C)P+ , Pf = [(1 - OV + rl]P + , Pi = (CV - rl)P 



+ 



Pi = -n = a/(0 2 - ml)(l - , (Pf) 2 = (if) 2 = Pf • (3) 



The polarization vectors e((/>) of the 4> meson are obtained from the orthogonality e(0) • Pi = and from the normal- 
ization e((j)) 2 = —1. The exact expressions are given, in the light-cone coordinates e ~ (e + , e~, t x , e y ), by 



1 / C[C+(l"0»? + i]-2^ [(1 - OV + f%][C+(l- Ov + rl] -2 



4\<t>) 



r t> \V2^[{+(l-()ri + rl] 2 -4rl v^lC + (1 - C)v + r 2 ] 2 ±r 2 
l(+(l-C)v + rl}^/(Ol-r 2 )(l-0 
^[C+(l-Ov + rD 2 -4r 2 



+ (1 - OV + rl? 4r 2 ' + (1 - C)V + r\\ 2 - ±r 2 

[C-d-Ov-r 2 ] 
^K+(l-Ov + rl] 2 -4r 2 

<%\<l>) = (0,0,0,1). (4) 

The terms proportional to will be neglected eventually. The kaon is treated as a massless particle. The <f> meson 
emitted from the weak vertex then carries the momentum P 3 = (ms/\/2)(0, 1 — 77, Ot)- Another equivalent, but more 
general, representation of e(0) is given by 

e T {4>) = , e T (*) = — ^— , (5) 

with p = Jp 2 + (P+ - Pf) 2 /2. 
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The three-body B meson decays are dominated by the contribution from the region, in which the <f>K pair possesses 
the invariant mass w 2 ~ 0(Am#) 0, A representing a hadronic scale. The orders of magnitude of the components, 



P+ ~ 0(m B ) , P~ ~ 0(A) , P T ~ OiVAms) , (6) 
are then implied. It is easy to obtain the power counting rules of the polarization vectors from Eq. Q , 

e+fo) ~ -^O(l) , effa) ~ -^O (jA/m B ) , e£(0) ~ -0(A/m B ) , 

eW + (0)~eW-(0)~of JIf^\ , e [ T )x {<t>) ~ 0(1) . (7) 



In the heavy-quark limit the hierarchy P + 3> Pr S> P _ corresponds to a collinear configuration, and suggests the 
employment of the new nonperturbative inputs, the <pK two-meson distribution amplitudes. For the PP system, there 
is only a single twist-2 distribution amplitude associated with the structure 7^, and two two-parton twist-3 distribution 
amplitudes associated with the structures / (the identity) and a^ v 0,0, E3- Here a higher- twist distribution amplitude 
means that its contribution is suppressed by powers of w/ tub - For the VP system, the relevant structures are more 
complicated: three twist-2 distribution amplitudes are associated with 7^75 and 07^75, and hve twist-3 distribution 
amplitudes with 7^75, 07,1,75, 75 and 7^. To decompose the two-meson distribution amplitudes into the components 
of different twists, we introduce the polarization vectors of the <pK system, 

e i = -L (1,-^,0,0), e£> = (0,0,1,0), e^ 2) = (0, 0, 0, 1) . (8) 

A two-pion distribution amplitude has been related to the pion distribution amplitude through a perturbative 
calculation of the process 77* — ► tt + tt~ at large invariant mass w 2 |3*cj . In this work we adopt a similar trick: we 
calculate perturbatively the matrix elements, 

( ( f>(P 1 ,e^))K+(P 2 )\u(y-)Ts(0)\0) , (9) 

using the (f> meson and kaon distribution amplitudes up to twist 3 |3ll l32j| , where T represents a structure among 7, 
75: Ifi > 7/x75 and 07^,75. The matrix elements can be expressed as the products of the corresponding form factors with 
the kinematic factors. For example, the matrix element for T — 7^75 is written as the product of the form factor Py 
with the kinematic factor (Pi — Pa)^- The kinematic factors are then approximated in terms of the momentum P and 
the polarization vectors e of the <j>K system according to the power counting rules in Eqs. © and 10. The resultant 
(^-dependent coefficients in the approximation contribute to the £ dependence of the <pK two-meson distribution 
amplitudes. 

We then derive the decomposition up to Oiw/ms), 

(<l>K+\u(y-)^ 5 s(0)\Q) = P M / dze izP -y$ {] (z,{,w) , (10) 







^K + \u( y -)a^ l5 s{0)\Q) 



-1 |(e TM P, - e Tv P„) J dze lzP - y <S> T (z, C, w) 

+-(Pl»P2 V - PluP2») [ dze izP - y <S> 3 (z,(,w)X , (11) 
w Jo ) 

(<PK+\u(y-h 5 s(0)\0) = w [ dze izp -^ p (z,C,w) , (12) 

Jo 

(^K+\u(y-)^ S (0)\0) = irJU—e^^pPnZ. [ dze lzP -y<$> v (z, £ w) , (13) 

f -n- J 

(^if+|«(l/-)Zs(0)|0) - 0, (14) 



where z is the momentum fraction carried by the spectator u quark, and n_ = (0,1, Ot) a null vector. We have 
adopted the convention e 0123 = 1 for the Levi-Civita tensor e^P s . The above decomposition applies to other VP 
systems, such as K*tt, pK, 

Below we present some details of the expansion of the kinematic factors. For Eq. i|10fl , we have applied 

(Pi-P 2 ) M «(2C-1)P M , (15) 
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where the coefficient 2<^ — 1 is absorbed into the distribution amplitude giving its £ dependence. Similarly, we 
have approximated the kinematic factor for the matrix element in Eq. I|f f fl . 

CTuMPi* - e Tv {4>)P^ « C(eo>P* ~ e T uP„) , (16) 

where the coefficient £ is absorbed into <&t, and is a transverse polarization vector of the <pK system. The 
contribution from another distribution amplitude $3 can be combined with that from <$>t via the approximation, 

I (Pi„P-Zu ~ PiuPiJ « 2v/C(W)(e^P v - e^,) , (17) 

where the coefficient ^/ £(1 — £) comes from Pf in the — > limit. Since the branching ratio is a sum over the 
transverse polarizations and e^, we omit the coefficient 2, and replace by the two possible ej^i- We have 
employed the approximation for the matrix element in Eq. (|f 3fl . 

2 wj 

-W^rMPfPf « 7T C^, P ^ T P p n a _ . (18) 



For this structure, the <p meson emitted from the weak vertex must carry a transverse polarization, and a non- vanishing 
hard kernel demands that the subscript \i denotes a transverse component. The coefficient £ is then the sum of £, 
C - 1, and 1 - ( from the combinations '(</>) = e^ 1)J "(0), P p = P+, P% = P 2 ~~), (e^W) = eJ' 1 ^, P p = P~, 
Pg = P})t and (e^ u ((f>) — e^p (</>), P p = P + , Pf = Pg"), respectively. A coefficient 2 for the last combination has 
been omitted for the same reason. 

Our strategy does not provide the z dependence. Assuming the z dependence of each C, w) to be asymptotic, 
we propose the parametrization, 



$n(«,Ci«0 


= — z(J — 
\/2lVc 


«)(2C"1) , 




3PtM ri 
= — z 1 - 


*)C, 




= 3i ^(l- 


*), 




= 3F >),(1- 
V2Ac 






= 3F ^z(l- 
V2Ae 


z)C. 



(19) 



The time-like form factors Pi 



lu) define the normalization of the <f)K two-meson distribution amplitudes. Note 



that these form factors are normalized to P||.t.3,p,u( to 0) = 1 in order to respect the kinematic threshold of decay 
spectra. Our strategy also reveals the power behaviors of the form factors in the asymptotic region with large w, 
Fu T (w) ~ 1/w 2 and F3 tPtV (w) ~ mo/w 3 , mo m 1.7 GeV m m being the chiral scale. Therefore, we further 
parameterize the form factors in the whole range of w for the evaluation of the nonresonant contribution: 



PiiM 

Ft(w) 
F 3 (w) 
F v (w) 



rut 



(w — m^) 2 + rrVu 



2 ■ 



(w - m^) 2 
F P H = 



,2 ' 



mom, 



(w - m^) 3 + m m| 



(w - m ) 3 + m m 2 T 



(20) 



where the two free parameters mil t, expected to be few GeV are determined by the fit to the measured r — > (j)Kv 
and B — * <f>K^ branching ratios |2q|. The form factors depending on the parameter my (my) are associated with the 
longitudinally (transversely) polarized <fr meson. 
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We stress that Eqs. (1101) and (|llfl contain not only the twist- 2 distribution amplitudes, but the twist-3 ones. The 
expansion in Eq. (|15f) corresponding to the component /i =_L generates 



we Tp [ dze lzP - v $ a (z,(,w) , 
Jo 

3F,|(w 2 ) 



/2Nr 



z(l-z)VC(W) 



(21) 



Similarly, we extract two twist-3 distribution amplitudes from Eqs. i|16|) and (|17|) corresponding to the components 
H,v = +,-, given by 



(e L „P v - e Lv P„) f dze lzP -y [(2£ - 1)$ 3 («, C, w) - $ t (z, C, 
Jo 



(22) 



For the </>-K" system, the above twist-3 distribution amplitudes lead to smaller contributions compared to <J> P and 
and have been ignored: because of m^/w 2 ~ 1, the range in Eq. (|28J) below indicates £ ~ 1, and that the contribution 
from $ a is suppressed by the factor ^/l — (. There exists a strong cancellation between $3 and <!> t in Eq. Q22[l. For 
other systems, such as pif, these twist-3 distribution amplitudes could be numerically important due to m 2 /w 2 <C 1 
in this case. 

For the B meson distribution amplitude, we use the model |l6[ . 



$b(x) = N B x 2 (l - x) 2 exp 



1 / XWlB 

2 l~^T 



(23) 



with the shape parameter = 0.40 ± 0.04 GeV j34[, and the normalization constant Nb being related to the 
decay constant Jb = 190 MeV (in the convention f„ = 130 MeV) via J Q <& B {x)dx = /s/(2y / 2^). The range of 
LOB is determined from a fit to the values of the B — > ir form factor from light-cone sum rules l35l |3(l . The above 
<f>s is identified as $ + of the two leading-twist B meson distribution amplitudes <I>± defined in [371 138|. Equation 
(|23|l . vanishing at x — > 0, is consistent with the behavior required by equations of motion [3^. It has been shown 
that the B meson distribution amplitude is normalizable in J. factorization theorem H, contrary to the conclusion 
drawn in the framework of collinear factorization theorem |4lL l42j . Another distribution amplitude $ b , identified as 
<3>b = ($_ — $ + )/v / 2 with a zero normalization, contributes at the next-to-leading power of A/ms Hi|- It has been 
verified numerically |43l | that the contribution to the B — > ir form factor from <&b is much larger than from &b- 
In summary, we calculate the hard kernels by contracting the quark-level diagrams with the matrix elements, 



(0\b(0)id(y-h\B(PB)) 

(<f ) K(P,e L )\u(y-) 3S (Q) l \0) 

(0jr(p,e T )|«(ir)j*(o)i|o) 



1 



/ 2Nc 
1 

1 



P-n. 



dxe- lxP - y [{/P B + m B h 5 ]ij$B(x) , 

dze izPy {( 75 f) lj $\\(z,C,w) + ('y B )ljW*pfa<,w)} » 

e^ P A^he u T P p n^ v (zX,w) 



(24) 



which follow Eqs. (|10() - (|13|1 . The calculation of hard kernels is as simple as of two-body decays. It is observed that the 
distribution amplitudes $||,r,3 give leading contributions, and those from $ Pjt , are suppressed by a power oiw/rriB- 
The r — > <f)Kv differential decay rate in the (f>K invariant mass is written as 



dT 

dw 



with m T being the r lepton mass, and V U8 the Cabibbo-Kobayashi-Maskawa matrix element. The B 
spectrum is written as 



dT 

dw 



r 2 ™ 4 



dC\M((,w)\ 2 , 



(25) 
(f>Kj decay 

(26) 
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with the amplitude, 



M(C,w) 
A(C,w) 



—V ts V tb m b A{(,w) , 

SttC Fm 2 B eT{^) ■ £t / dxidz — 
Jo %i 



P4 



(1 + z) ($ T (z, C, w) + $ 3 (z, C, w) VC(1 - 

(4 2) )c ? t(t( 2 ))l 



v ^(l-2z)$ t ,( z ,C,H 



zm ?, + P£ 



x\m 2 B 



(27) 



e(7) and Qj, represent the photon polarization vectors and the photon momentum, respectively, mj is the b quark 
mass, and Cf^ the corresponding effective Wilson coefficient [44|. All the terms of O(ry) in M. have been neglected 
for consistency. The requirement P T > leads to the bounds of £ as shown in Eq. l(2T)|) , 



ml/w 2 < C < 1 . 

The hard scales are chosen as the maximal virtuality in each quark-level diagram 0, IT^ | , 



P| 



t (2) 



x\m 2 B 



P 2 T 



(28) 



(29) 



The above collinear factorization formula is well-defined, since the invariant mass of the two-pion system, appearing 
through Pt, smears the end-point singularities from z — > 0. Even if one adopts a model of the B meson distribution 
amplitude, which vanishes only linearly in x\, Eq. (|27J) is still well-defined due to the presence of Pt- 

Because there exists only an upper bound for the measured r — > <j)Kv branching ratio, we also consider the 
r — > K*nv branching ratio, when constraining the parameter win. That is, we assume that mil in the two decay 
modes, i.e., the time-like <f>K and K*n form factors, do not differ much. The experimental data and the theoretical 
prediction from chiral perturbation theory are 



B(t 4>Kv) < 6.7 x 1(T 5 
B(t -> K*nv) ~ ° 7 " 1ft " 4 
B(B -> cj)Kj) 



26 

2.7 x 10" 4 lEl , 
(3.4 ±0.9 ±0.4) x 10~ 6 



(30) 



For the application of our formalism to the r decays, we shall trust it up to the order-of-magnitude accuracy, since 
the ratio r\ = w 2 /m 2 in this case is not a small parameter. With m\\T sa 2 GeV, we obtain from Eqs. f25J) and (|27|l . 



B(t -> 4>Kv) = 6.8 x 10~ 5 , 
B(t -> iTvn/) = 4.5 x 10~ 4 , 
S(S - <^T 7 ) = (2.9±g:J) x 10- 6 



(31) 



consistent with Eq. I|30[l (up to order of magnitude for the r decay data as stated above) . Our results are stable with 
respect to the variation of mil and m<r around few GeV. Therefore, the theoretical error in Eq. (|31|) comes from the 
variation of the shape parameter wg, which can be regarded as an estimate of the uncertainty from hadronic dynamics. 
The predicted B — > (frK'y decay spectrum is shown in Fig. ^ which exhibits a maximum at the <pK invariant mass 
around 1.3 GeV, consistent with our power counting rules. 

After constraining the two-meson distribution amplitudes, we predict the B — > (fxfrK decay spectrum in the <pK 
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FIG. 1: B — * cj>K"f and B — ► (jxj>K decay spectra in the (j>K invariant mass. 



invariant mass. For this mode, the amplitude M. is written as 



M = 



P(s) 



i=3 



+ e T ■ e ZT {6r)T^ s) 

$b(xi)$\\(z,C,w) 
X\zm B + Pj, 

x { [(1 + z)$|| (z, C, w) + 0?(1 ~ 2z)%(z, C, w) 



(32) 



8itCf17i b / dx\dz 



a s (t 



(*W) 



ZTTii 



P| 



P( S ) 



-2 v ^$p(z,C,' l 



8r<pTrCFTn B 



, g s (^K (s) (4 2) ; 

' 2 

dx\dz- T ""' 1 



}■ 



(33) 



x\zm B + Pf 



x | |$ x (z, C, ttf) + $ 3 («, Cj w)VC(l-C) + Zy/rj9 v (z, C, w) 

IlITln J 



a,(tl 1) )ai ,) (t' 1 ') 

Zfflj + Py 



(34) 



where £3t(</0 denote the polarization vectors of the 6 meson emitted from the weak vertex. The definitions of the 

Wilson coefficients a\ (i) are referred to [47^. For a similar reason, we have dropped all the 0(rf) terms. Equations 
(|27|l , and (|34|l represent the amplitudes of the B meson transition into a VP meson pair associated with different 
effective operators. We display the predicted B^ — > cjxpK decay spectrum in Fig. ^ which also exhibits a maximum 
at the 6K invariant mass around 1.3 GeV. Integrating the spectrum over 77, we obtain the branching ratio without 
the resonant contribution in the <bd> channel, 



B(B ± -> 66K*) = (1.3l£S) x 10 



(35) 



The uncertainty arises from the variation of the shape parameter lob of the B meson distribution amplitude. 

We have examined other sources of theoretical uncertainty. The correction to the branching ratios from the ne- 
glected 0(77) terms is about 10%. To investigate the uncertainty from different parametrization of meson distribution 
amplitudes, we have tried 



<&' B (x) = N' B x(l - x) exp 



1 / xrriB 
2 



(36) 
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First, the shape parameter uj' b = 0.9 GeV is determined from the fit to the value of the B — > tt transition form factor 
about 0.3. The model &' B (x) is then employed to fix the <j)K two-meson distribution amplitudes from the data of the 
B — > <f>Kj branching ratios. It is observed that the symmetric z dependence in Eq. (|19|1 should be modified into 

2(1 -z)->z(l-2)[l + 0.5(1 - 2z)} , (37) 

which is reasonable since the <j> meson is heavier than the kaon. After going through the above procedure, we predict 
the B — > <p4>K branching ratio using the distribution amplitudes in Eqs. (|36|) and l|37|) . and find that the result 
increases only by 8%. We have also checked the sensitivity of our prediction to the parametrization of the time-like 
form factors. Obeying the normalization and the asymptotic behavior required by PQCD, the models with (w — m^) 2 
[(w — to^) 3 ] being replaced by w 2 —rr& [(w 2 — m|) 3 / 2 ] are also allowed. Adopting the B meson distribution amplitude 
in Eq. (|23[) . the t — > <j)Kv and B — * (f)K~/ data just imply a slight increase of the parameters m» and mj- to 3-4 
GeV. Then we predict the B — > <fxf>K branching ratio using the new parametrization, which is enhanced only by 12%. 
The above investigations indicate that the PQCD predictions will be insensitive to the parametrization of meson 
distribution amplitudes, if the procedure of determining meson distribution amplitudes is followed. 

Note that the B ± -> 00X ± branchin g ra tio has been measured to be B{B ± -> 00X ± ) = (2.6±o g ± 0.3) x 10~ 6 
for a 4>(j) invariant mass below 2.85 GeV [27]]. We suggest that the decay spectrum in the cf>K invariant mass should 
also be measured (only the spectrum in the <jxf) invariant mass was presented in 27] ) , such that the dynamics of the 
B — > VP transition can be explored. To derive the spectrum in the <j><j> invariant mass, we need to define the VV 
two-meson distribution amplitudes, which will be discussed in the future. 

We thank H.Y. Cheng and Y.Y. Charng for useful discussions. This work was supported in part by the National 
Science Council of R.O.C. under Grant No. NSC-92-2112-M-001-003 and No. NSC-92-2112-M-006-026. 
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